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Abstract: This note provides an alternative approach to the momentum decay and 
thermal evolution of decoupled massive particles. Although the ingredients in our 
results have been addressed in Ref. U, the strategies employed here are simpler, and 
the results obtained here are more general. 
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1. Introduction 

As is well known, for the freely traveling massless particle like photon in an expanding 
FLRW universe, the frequency or energy will vary inversely proportional to the scale 
factor, which implies that the number density of massless particles still keeps its thermal 
spectrum form with a redshifted effective temperature although these particles went 
out of the thermal equilibrium into the free expansion as time passed. This is the 
physical foundation for the cosmic microwave radiation background currently observed 
by us. Now a natural question arises, namely, does the above fact also apply to the 
massive particle? Not only does this question possess a theoretical interest by itself, 
but also acquires a practical implication in cosmology since neutrinos and antineutrinos 
are believed to be massive. However, to my best knowledge, this issue has not been 
addressed in literatures except in Weinberg's cosmology book published recently]]]]. 

The purpose of this note is to provide an alternative approach to this issue. The 
strategies employed here are simpler, but the results obtained here are more general. 
Notations and conventions follow Ref.pJ. 

2. Momentum Decay 

In general curved spacetime, a particle of mass m freely travels along the timelike 
geodesic t]{t) with r the proper time, which means that U a = (J^) a gives the geodesic 
equation 



U a V a U' 



b 







(2.1) 
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with U a U a = — 1. Assume there to be a family of observers Z a along the geodesic, then 
we have 

d 4- = U a V a (-mU b Z b ) = -mU a U b V a Z b , (2.2) 

GST 

where E = —mU b Z b is the energy of massive particle measured by the observers. 
Now for the expanding FLRW metric 

dv 2 

ds 2 = -dt 2 + a 2 {t) i 1 _ Kr2 + r 2 {dd 2 + sin 2 Odip 2 )] (2.3) 

with K = 1,0,-1 for closed, flat, and open universes respectively, if the observers are 
chosen to be the isotropic ones as usual, i.e., Z a = (J^) a , we have 

V a Z b = -h ab , (2.4) 
a 

where the dot denotes the derivative with respect to the time t, h ab is the induced 
metric on the surface of constant t, given by h ab = g ab + (dt) a (dt) b . Plugging Eq.( [2.4| ) 
into Eq. (|2~^ ), we obtain 

= -m-U a U b h ab = -m- -1 + (U a Z a ) 2 } = , (2.5) 

clt a a ma 



which implies 



da_ m dt dE _ EdE _ 1 d(E 2 - m 2 ) _ dp 
a dT E 2 — m 2 E 2 — m 2 2 E 2 — m 2 p ' 



where p = \JE 2 — m 2 is the magnitude of momentum of massive particle measured by 
the isotropic observers. Whence we know that for a freely traveling massive particle in 
an expanding FLRW universe, it is its momentum rather than energy that goes like 1 

poci. (2.7) 
a 

It is noteworthy that this result is also obtained in Ref.Q, where, however, the method 
employed seems somewhat complicated, and some approximations are also made. 



3. Thermal Evolution 

Let us assume that during the evolution of our universe, there exists a last scattering 
surface at the time ti when some kinds of massive particles such as neutrinos and an- 
tineutrinos suddenly went from being in thermal equilibrium to a decoupled expansion. 

1 Of course, the momentum of a massless particle shares the same behavior since its momentum 
equals energy. 
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Then according to Eq. (|2.7|) the massive particle that has momentum p at a later time 
t would have had momentum p L = p-^^j at the time t^. So the number density of 
massive particles at the time t with momentum between p and p + dp would be 

n(p,t)dp = (^^-) 3 n(p L ,t L )d(p L ) 

a(t L ) 3 4Trgp 2 L dp L 1 



a(t) (27rh) 3 exp [(v/pl + m 2 - fi d )/kT d ] ± 1 
4:7rgp 2 dp 1 

(2tt£) 3 exp [{y/pi + m 2 - fi d )/kT d ] ± 1 
4ngp 2 dp 1 

(2vrfr) 3 exp [(^p 2 + m\ - fi e )/kT e ] ± 1 ' 



(3-1) 



Here the factor (^^-) 3 in the first step arises from the dilution of particles due to the 
cosmic expansion. The Fermi-Dirac and Bose-Einstein distributions are employed in 
the second step, where g is the number of spin states of the particle and antiparticles, 
\i d and T d denote the chemical potential and temperature in thermal equilibrium at the 
last scattering surface, respectively, and the sigh is + for fermions and — for bosons. 
We introduce the effective mass, chemical potential, and temperature in the last step, 
i.e., m e = m^lp // e = n d ^^-,T e = T d ^jy. Therefore the form of the Fermi-Dirac 
and Bose-Einstein distributions are preserved for the thermal evolution of decoupled 
massive particle, with the effective mass, chemical potential, and temperature varying 
inversely proportional to the scale factor a at the same time, which implies that the 
ratios among the effective mass, chemical potential, and temperature remain constant, 
just as before decoupling. 

Note that just by taking the mass to be zero the above argument obviously reduces 
to the massless case, where the result is also obtained in Ref. |l| by the thermodynamic 
method rather than the simpler dynamical picture employed here. In addition, although 
the spectrum has still kept the form of the Fermi-Dirac and Bose-Einstein distributions 
since decoupling, it is not the thermal spectrum with the effective temperature and 
chemical potential since the effective mass is not equal to the static mass. The unique 
exception is the massless case. 
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